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Intensive high-pressure experiments on 
the transition-metal monoxides have re- 
vealed that FeO fl], ||, |, g § and MnO |§, 
0, ^ undergo pressure-induced phase tran- 
sition. The high-pressure phase of FeO 
was identified as the NiAs (B8) type [| 
while that of MnO is yet unclear 
The present theoretical study predicts 
that the high-pressure phase of MnO is a 
metallic normal B8 structure (nB8), while 
that of FeO should take the inverse B8 
structure (iB8). The novel feature of 
the unique high-pressure phase of stoi- 
chiometric FeO is that the system should 
be a band insulator in the ordered an- 
tiferromagnetic (AF) state and that the 
existence of a band gap leads to spe- 
cial stability of the phase. The observed 
metallicity |l^] of the high-pressure and 
high-temperature phase of FeO may be 
caused by the loss of AF order and also 
by the itinerant carriers created by non- 
stoichiometry. Analysis of x-ray diffrac- 
tion experiments provides a further sup- 
port to the present theoretical prediction 
for both FeO and MnO. Strong stability of 
the high-pressure phase of FeO will imply 
possible important roles in Earth's core. 

FeO and MnO at normal pressure are typical 



examples of the Mott insulator [ 1 1 1 and their ba- 
sic properties are governed by the electron corre- 
lation whose strength is measured by U /W with 
U the effective Coulomb interaction integral be- 
tween d electrons and W the d band width. As 
the d band width is most directly controlled by 
pressure, the high-pressure studies have been re- 



garded clS cl useful way to understand the basic 
properties of Mott insulators. Besides, the high 
pressure phase of FeO may have an important 
implication in the earth science because FeO is 
believed to be one of the major constituents of 
earth's deep mantle. 

Under normal pressure and room tempera- 
ture, FeO and MnO both take the rock-salt (Bl) 
structure with rhombohedral distortion along 
the < 111 > direction. This distorted Bl struc- 
ture is called rBl hereafter. It was observed that 
FeO undergoes a phase transition from the rBl 
structure to a high pressure phase at pressures 
above 70GPa [Q, ^, whose crystal structure was 
recently assigned to be the B8 by the analysis 
of x-ray diffraction peak positions |Q, |5|. On the 
analogy of almost all of the compounds with the 
B8 structure, a natural idea for the B8 (NiAs) 
FeO may be such that Fe occupies the Ni site 
and O the As site. This structure is named nB8. 
However, another structure, which is named iB8, 
is possible by exchanging the Fe and O posi- 
tions [^]. More recently, not only shock com- 
pression but also static compression |^] ex- 
periments on MnO also showed the existence of a 
possible high-pressure phase above 90GPa. How- 
ever, the crystal structure has not been deter- 
mined yet. 

In the present work, we have performed the 
first-principles density functional calculations to 
give valuable hints to such problems as those 
mentioned above. The electron-electron inter- 
action is treated by the generalized gradient 
approximation (GGA) |^] like in other simi- 
lar calculations [|l^, However, the GGA is 
not powerful enough to treat the strong corre- 



lation system such as Mott insulators. There- 
fore, the GGA calculation is supplemented by 
the LDA+U method [15| particularly in the low 
pressure regime where MnO and FeO are re- 
garded as Mott insulators. 

Highly converged total energy calculations are 
performed for different crystal structures (Bl, 
B2 (CsCl type), nB8, iB8) and spin structures 
(ferromagnetic (FM), AF) for different volume 
(shown in Figure 1). Surprisingly, the iB8 struc- 
ture with the AF ordering is the most stable for 
FeO among several structures (including rBl) in 
the whole volume range. First we pay atten- 
tion to the relative stability between the IBS and 
nB8 structures in the compressed volume range, 
where the GGA calculation will be reliable be- 
cause of the reduced U /W. Clearly for FeO, the 
IBS structure is more stable than that of nB8. At 
this stage, two fundamental questions are to be 
answered: 1) why is the nB8 structure is realized 
rather than iB8 for most of the transition-metal 
compounds with the B8 structure (including the 
high-pressure phase of MnO, see Fig. 1(b))?; 2) 
what is special for FeO for the strong stability of 
iB8? Anions around a transition metal ion form 
an octahedron in nB8 and a trigonal prism in 
iB8. In the latter case, absence of the local in- 
version symmetry about the transition metal ion 
reduces the strength of the hybridization of the 
3d orbitals with the oxygen 2p orbitals. As the p- 
d hybridization contributes to the stability of the 
structure, this aspect favors the nB8 structure 
and will explain the general feature that the nB8 
structure is actually realized in most cases. How- 
ever, the calculated energy level diagram shows 
existence of a well defined band gap between the 
conduction band and the valence band for the 
AF-iB8 structure of FeO in the whole volume 
range shown in Fig. 1(a), contributing to special 
stability of this structure. 

A puzzling feature of Fig. 1(a) is that the iB8 
structure is significantly more stable than the 
rBl structure even at equillibrium volume at zero 
pressure. This reflects the fact that the GGA 
cannot describe the electronic structure of Mott 
insulators properly: the GGA incorrectly makes 
FeO at normal pressure metallic. In order to 
reproduce the correct ground state of FeO, we 
have to take account of the local electron cor- 
relation and the spin-orbit interaction as well. 



We adopt the LDA+U method with the LMTO 



(Linear- Muffin-Tin Orbital) basis |16] to perform 
calculations with these ingredients. A reasonable 
value of 4 eV as the effective Coulomb interaction 
parameter U^s makes the rBl structure insulat- 
ing and more stable than the iB8 structure at 
normal pressure. On the other hand, for a com- 
pressed volume with experimentally determined 
lattice parameters [||, the same calculation still 
makes the iB8 structure most stable. Therefore, 
the puzzling feature of Fig. 1(a) is removed by 
the electron correlation effect. 

According to the present calculation, the high- 
pressure phase of FeO with the iB8 structure 
should be insulating in contrast to the experi- 
mental observation of metallic behavior j^, 10 1. 



There are two possible origins in this disagree- 
ment. Our calculation assumes the stoichiomet- 
ric FeO, while real samples contain about 5 % Fe 
deficiency. As the iB8 FeO is a band insulator 
rather than a Mott insulator, itinerant carriers 
will be doped by Fe deficiency. Another possi- 
bility is related to the magnetic disorder in the 
temperature range of the observed metallic be- 
havior (~ 1000 K). According to the present cal- 
culation, even if the crystal structure is the iB8 
type, the FM order makes the system metallic, 
suggesting that the metallic state may be caused 
by the loss of the AF order at such a high tem- 
perature range. 

There is another strong evidence for the iB8 
structure as the high-pressure phase of FeO. We 
have found that the intensity profile of the ob- 
served x-ray diffraction pattern ||5| can be re- 
produced only by the iB8 type but not by the 
nB8 type. To be more concrete, we consider 
the relative intensity between (100) and (101) 
peaks as an example. Experimentally, the lat- 
ter is stronger than the former. This feature is 
correctly reproduced only by the iB8 structure. 

As for MnO, the present calculation shown 
in Fig. 1(b) suggests that the most stable high- 
pressure phase will be the nB8 structure rather 
than the B2 Q and iB8 structures. Detailed 
comparison of the total energies predicts that the 
FM-nB8 structure has the lowest energy rather 
than the AF-nB8 structure, though the energy 
difference is rather marginal. 

The present theoretical prediction can explain 
the recent experiment S on MnO in a high pres- 
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Table 1: Observed [8] and fitted x-ray diffraction 
patter of MnO at 137GPa after laser annealing. 
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t dcxp and dfit are experimental and fitted d- 
spacings respectively in unit of A. The fitted 
hexagonal unit cell has a = 2. 441 A and c = 
5.076A. 

I The relative intensities of the peaks are de- 
scribed as strong(s), medium(m) and weak(w). 

sure range (> 120 GPa) very well. The con- 
sistency of the assignment of the experimental 
x-ray diffraction peaks by the nB8 structure is 
demonstrated in Table 1. Almost an exact fit of 
peak position and good agreement of the inten- 
sity profile can be obtained except one peak for 
dcxp = 1.844 A which may originate from some 
intermediate phase. The intensity of this peak 
is actually reduced after annealing. The volume 
(=7.9 cm^/mole) and the c/a{= 2.08) estimated 
by fitting the peak positions are in good agree- 
ment with the present calculated results. 

Figure 2 shows the P-V relation for MnO ob- 
tained by the present calculation and also by ex- 
periments. Clearly, in the low pressure range, 
the sample is in the AF-rBl phase. The ex- 
perimental data point around 120 GPa by shock 
compression is located just in between the curves 
corresponding to the rBl and nB8 structures, 
suggesting that the sample is in a mixed phase. 
On the other hand, the highest pressure data 
point obtained by the static compression fol- 
lowed by laser annealing is just on the line of 
the P-V curve for the nB8 structure. 

We would like to make a brief comment on 
the c/a value of the nB8 of MnO and iB8 of 
FeO. In the high-pressure phases, c/a of both 
systems exceeds 2.0, being unusually large com- 
pared with the values for most of other related 



systems. However, we have found that c/a is 
an increasing function of Tc/ra with ra (re) de- 
noting the anion (cation) ionic radius and that 
the c/a values of the MnO and FeO are on the 
extrapolated line of this general trend. 

Finally, as one of the important consequences 
of the present analysis, strong stability of the iB8 
phase of FeO may imply strong tendency of the 
incorporation of oxygen in Earth's core. 

Methods 

Pseudopotential calculation. The plane- 
wave basis pseudopotential method is used to 
perform the structural optimization efficiently. 
The 2p states of oxygen and 3d states of Mn 
and Fe are treated by the Vanderbilt ultrasoft 
pseudopotential |l^. The cutoff energy for the 
plane-wave expansion is 4.9 x 10^ eV, and sam- 
pling of k-points gives absolute total energy con- 
vergence to better than 10^^ eV/fomular unit 
(~ 1.0 KJ/mole). Moreover, note that the con- 
vergence in the relative energy between differ- 
ent structures is generally order of magnitude 
much better than the convergence in the abso- 
lute total energy. The equation-of-state param- 
eters obtained by the present calculations agree 
well with those obtained by other similar calcu- 



lations 1 14, 15]. High credibility of the present 
approach has been well confirmed. 

LDA-|-U method. The interaction between 
localized d electrons is explicitly taken into ac- 
count through the screened effective Coulomb in- 
teraction parameter C/eff- 

Estimation of the p-d hybridization. The 
reduction in the p-d hybridization in the iB8 
structure compared with that in the nB8 struc- 
ture was confirmed in the present calculation 
by estimating the energy separation between the 
centers of the p and d bands. 
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Figure 1: The total energies of FeO (a) and MnO (b) as a function of the volume. The least- 
squares-fitted curves to the Murnaghan's equation of state [19] are shown. 
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Figure 2: The Pressure- Volume curves for different phases of MnO. 
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